Degenerate homogeneous parabolic equations 
associated with the infinity- Laplacian 



Manuel Portilheiro* 
Juan Luis Vazquez''' 

Abstract 

We prove existence and uniqueness of viscosity solutions to the degenerate parabolic 
problem ut = AJJ^m where Ajj^ is the /i- homogeneous operator associated with the infinity- 
Laplacian, Ajj^it = \Du\^~^ {D^uDu, Du) . We also derive the asymptotic behaviour of u 
for the problem posed in the whole space and for the Dirichlet problem with zero boundary 
conditions. 

1 Introduction 

In this paper we study the following class of degenerate parabolic equations for u : Q — !■ R, 
u = u{x, t) with a parameter h > \: 

(1.1) ut~l^l^u = 0, inQ, 

(1.2) u = g, onr, 

where denotes the ^.-homogeneous degenerate operator, 

d 

A^u := \Du\^-^ u^,x,u^,u^y 

Note the h stands for the homogeneity degree of the operator. The equation is posed for (x, t) S 
Q := X (0,r); where il C M'' is open (not necessarily bounded), < T ^ oo, F = dpQ the 
parabolic boundary of Q and g : F — > R is a continuous function. 

Equation (|l.ip has been studied in two particular cases, h = \ and /i = 3, which correspond to 
the two common definitions of the infinity-Laplacian. For = 1 we find the usual 1-homogencous 
infinity-Laplacian operator, Aoo w = l-^^I^^Si 

that for other h we get the for- 
mula A'^u = |Z?u|''~^(Aoo u). The infinity-Laplacian operator has received a lot of attention 
in the last decade, notably due to its application to image processing. Most of the interest- 
ing properties of the cxD-Laplace equation can be found in Crandall's enjoyable paper [C] (the 
long list of references therein covers most of what we avoid listing here). Regarding the evolu- 
tion of the infinity-Laplacian, the above two cases {h = 1, 3) have been studied in a number of 
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works. The 1-homogeneous version has been considered by Wu, Juutinen-Kawhol and Barron- 
Evans- Jensen ( |W[ IJK| IBEJj : see jJKj for additional references in apphcations), whereas the 3- 
homogeneous equation has been treated by Crandall-Wang, Akagi-Suzuki, Akagi- Juutinen-Kajikiya 
and Larcngot-Stinner ( [CWl RSl [QKI IIB] ) . 

We think that equation (|l.ip . and its generahzation (jl.Sp . apart from their apphcations, have 
intrinsic interest, as they are degenerate parabohc, singular when /i < 3, do not correspond to a 
geometric equation, are not variational and are not in divergence form (with the exception of the 
1-homogeneous infinity- heat equation in two space dimensions, |DGIMR] ). They constitute a class 
of equations with particular difficulties and properties. 

It is also worthwhile pointing out that radial solutions of p.ip (or if we consider the equation in 
one space dimension) correspond to solutions of the p-Laplace evolution equation 

ut = div {\Du\P-'^Du) 

in one space dimension with h = p ~ 1. In fact, a comparison between the radial solutions in 
Section [3] and solutions of the p-Laplace evolution make this evident. 

Here we consider all cases h > 1. Since the results and techniques for /i = 1 are slightly different, 
we will not consider this case below and will treat it in a separate paper. For the case /i = 3 we 
extend slightly both the existence and the asymptotic results. Our existence result allows for an 
extra source term H{u) (see (|1.3|) bellow) and we do not impose the exterior sphere condition. 

We are specially interested in the asymptotic behaviour of solutions. Wc make the assumption of 
nonnegativity on data and solutions. For bounded domains we do not impose any extra conditions 
on fi. We also study the Cauchy problem posed in the whole K'', where the decay rate is different, 
as expected, and we get a precise convergence result to a self-similar solution of Barenblatt type. 



1.1 Main results 



To obtain the asymptotic behaviour on bounded domains it will be useful to prove existence for a 
slightly more general equation, allowing for zero-order terms. Therefore, we consider 

(1.3) Ut^A^u + Hiu), 

where : M — )■ ffi is a continuous function satisfying H{0) = and growing at most linearly at 
infinity, that is, there exists M > such that \H{y)\ ^ 

Theorem 1.1 (Existence). Let Q C M'*+^, g : F — > M 6e continuous and H : M. ^ M. be as above. 
If either Q is bounded or g is bounded, uniformly continuous in F and 

lim sup \g{x, t) \ — Q, 
R^°° |2-|5si?.,te[o,T) 

then there exists a unique continuous viscosity solution u of (|1.3p with u ^ g on F whose modulus 
of continuity depends on the modulus of continuity of g. If g is Lipschitz continuous, then so is u 
and its Lipschitz constant depends on the Lipschitz constant of g. 

Remark. In the case h = 2> this slightly extends jAS| Theorem 2.5], not only regarding the extra 
source term H(u), but as we do not impose the exterior sphere condition on the boundary of f2 
and we consider unbounded sets. 

We want to analyse the large time behaviour of the solutions of (|l.ip - (|1.2[) in two cases. The first 
is the Cauchy problem, Vl = W^, with compactly supported initial condition g = uq ^ Cc(M'*). 



(1.4) 



Ut = A'^u, in M'* X (0,oo), 

u{x, 0) = uo{x). 
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In this case u converges to a radial decaying solution. 

Theorem 1.2 (Asymptotic behaviour for the Cauchy problem). Let u he the unique viscosity 
solution of (|1.4p with uq S Cc(K'*) not identically zero, Up ^ 0. Then, for some i?* > 0, 



h+l _ h + 1 h+1 

R 1^ -t \x\ 1^ 



h-l 
+ 



lim t 2h sup \u{x,t) - Bfi hix,t)\ = 0. 
where B^ h *s the Barenblatt solution, which takes the form 

BR,h{x,t) := Qit^A 
as given in p.ip . 

Next we analyse the Dirichlct problem with homogeneous boundary conditions, 

fut-A^u^O, iuQ, 
(1.5) Iu{x,0)^uq{x), on VI, 

[ii(a;,t) = 0, for .T e t > 0, 

where uq e C{Vl) satisfies uq{x) = for x e dil, VL bounded and Q = x (0, oo). For nonnegative 
initial data, we obtain the following. 

Theorem 1.3 (Asymptotic behaviour for the homogeneous Dirichlet problem). Suppose Uq £ 
C'(fl), u ^ 0, uq ^ and uo{x) = for x G dfl. There exists a continuous function : — ^ R 
such that the solution of (jl.Sp . u satisfies 



lim t'^-iu{x,t) = Fn{x). 

Remark (1). Naturally, the result holds for constant boundary data with initial condition above 
(below) this constant, since adding a constant to u still gives a solution of the problem. 

Remark (2). We will show (Theorem 16. 3p that the function Gn {h — l)'^^Fn is a positive 
solution of the eigenvalue problem 

-ASoG = G. 

or 

-AooG ^ 



This is a variation of the eigenvalue problem for the infinity-Laplacian studied in [j] , 

- A^ G - AG 



and a variation of the nonlinear eigenvalue problem considered in |PV] 

-AooG = AGP, 0<p<l 



(in connection with this, [JR| considers "large solutions" of this equation with p > 1). Note also 
that the function {x,t) i— >■ t'^-'^ Fq_{x), which we call the friendly giant, is a generalization of the 
solutions of the form (13.31) for the radial case. 



Notations 

We always assume C is open and denote by Q a solid cylinder of the form Q x (io,ti), 
to <ti. We will also separate the parabolic boundary of Q into a lateral boundary and the base, 
dpQ = diQ U dbQ, with diQ = dVl x (to, ti) and d^Q = Q x {to}. 



3 



Contents 



1 Introduction [l| 

1.1 Main results [3 

2 Viscosity solutions Q 

3 Special solutions 

3.1 Similarity solutions 

3.2 Separation of variables 

3.3 Traveling waves 

4 Existence [l^ 

4.1 Lipschitz estimate in time [13 



4.2 Holder continuity in space [11 



4.3 Lipschitz estimate in space (e = 0) [13 

5 Asymptotic behaviour in the whole space [l^ 

5.1 Decay rate [l^ 

5.2 Asymptotic behaviour. Proof of Theorem ll.2l M 



6 Asymptotic behaviour on bounded domains [17 



6.1 Decay rate [17 

6.2 Friendly giant in n. Proof of Theorem [Ol [iS 

7 Comments and extensions [l9 



2 Viscosity solutions 

We start our analysis by defining viscosity solutions for Equation (|1.1|) . Note that since the equation 
is singular for h € (1,3) at points where the gradient of the function vanishes, the usual definition 
of viscosity solution needs to be adapted at the points of singularity. For the singular case we 
could use the definition from |CGG| (see also [JKj for the case h = 1). However, the singularity 
is removable, and since we are not considering here h = 1 we can simply recast Equation (|l.ip as 
follows: 

We can write this equation as 

(2.1) ut + F,,{D^u, Du) = with Fh{M,p) -\p\''-^{Mp) -p 

and Fh ■ SdX- (K'^ \ {0}) where Sd is the set oi dx d real symmetric matrices. Since, even for 
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1 < h < 3, lini|p|_i.o Fh{M,p) = for every M G Sd, we can take the continuous extension of Fh, 

Fh{M,p) 



Fh{M,p) if p 7^0, 
if P = 0, 



and the usual theory of viscosity solutions can be applied to F^- For simplicity we abuse notation 
and write F^ to denote its continuous extension F^- 

Definition 2.1 (Viscosity Solutions). An upper [lower] semicontinuous function u : Q — > R is a 
viscosity subsolution [supersolution] of (jl.ip in Q if and only if for every Pq £ Q and every 
function <f) G C^'^(Q) touching u from above [below] at Pq, that is, 

u{Po) = 0(Po) and 
■ u(P) < 0(P) [w(P) > 0(P)] for every PeQ\ {Pq}, 

we have 

(2.3) cj)t+Fh{D^(b,Dcj)) ^0 [4>t+Fh{D^4>.D4,)-;^(S\. 

A continuous function u G C{Q) is a viscosity solution when it satisfies both inequalities. 

Let us also record here the comparison property |CIL1 Theorem 8.2] that directly applies to 
([Illll-dlH) (see also [Ml Theorem 2.3]). 

Theorem 2.2 (Comparison Principle). Let u be an upper semicontinuous viscosity subsolution, 
and V be a lower semicontinuous viscosity supersolution of such that u ^ v on dpQ. Then 

u ^ V in the whole cylinder Q. 

Furthermore, if u and v are viscosity solutions of (|l.ip , then 

max \u ^ v\ ^ max \u — v\. 

Q dpQ 

Remark. Equation (jl.ip has a few symmetries which we will explore. Whenever u{x, t) is a classical 
solution, it is immediate to check that for any orthogonal dxd matrix O, for any (cCq, to) in R'' X R, 
c G R and any A and s positive, 

v[x, t) =u {0{x — Xo),t — to) + c 

(2.4) u)(a;,t) =Au(sx,A''"\s''+4) 

z(x, t) = — u{x, t) 

are also solutions in an appropriate domain. As usual, the argument can be transposed to viscosity 
solutions. 



3 Special solutions 
3.1 Similarity solutions 

We need to obtain certain explicit solutions. Let us first look at similarity solutions of (|l.ip . We 
are mainly interested in radial solutions, u{x,t) ~ v{r,t). It is easy to see that v should solve 
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which is a 1-dimensional p-Laplacian equation with p ^ h+\. Similarity solutions for this equation 
are well known. In fact, if we take v of the form 

v{r,t) = with y = rt^ . 

Then, we get for / the equation 

P-I [af{y) + (3yf{y)] = r'^+Z^C'+D {f{y)f~' f'{y). 
Choosing a = /? = — we obtain 

ayfiy) = I {ny)t + C. 
If we take C = we can integrate once more to get (where C is a diffcrcnc constant) 

I{y) 



C-['-y^-^y'^ 



which leads to the usual similarity solution 

(3.1) BR,h{x4) ^ Bh{x,t) := Cht-^ 



h+l h+1 , h + 1 



h-l 
+ 



where 



Ch 



2 [h + 1. 



Note that R denotes the radius of the positivity set of Bu h at time t = 1. It is straightforward to 
prove this defines a viscosity solution. 

Proposition 3.1. The function Bji^h defined in (|3.ip is a viscosity solution of p.ip in x (0, c»). 



Proof. Let us for convenience define A(r, t) = 
dtBh{x,t) : 



h+1 _ h+1 h.+ l 

R h ^ t 2h^r h 



We have 



Ch 2h + l , h Ch(h+1) 2h(h+l) + l h+l 1 

—t — —A~ + -P --t ^ r—A~, 

2h 2{h-l)h 



„„ , N Cft(/l + l) 2h+l 1 . 1 X 

DBh(x,t) = — ) ' t^^rT^A— — 
^ ' ' h-l \x\ 



2.x X 



D^B.ix,t) = j;^t~^rV^^A-^^ 

ciAh + l) 2h+i h-i . 1 fh — \x®x 
+ -^^T r-^t —A~ I — : — / 



h~l 



h 



From this computation we see that the set of points where B^ is not twice differentiable is 
W := {(a;,i) £ R"^ x (0, oo) | ^(|x|,i) = or a; = O}. Then, from the construction, B^ is a clas- 
sical solution of (|l.ip in R'' x (0, oo) \ W . To check B^ is a viscosity solution we just have to verify 
at the points of W . Note also that Bh is everywhere and DBh = on W . 



Let Pq = (s^Oi^o) £ be arbitrary and assume (j) touches u from above at Pqi that is, it satisfies 
(1221). Since Bh is C\ it must be that 0t(-Po) = dtBh{Po) and D(t){Po) = 0. If xo ^ 0, then 
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_2h±l h+1 

^{\xo\,tn) = and (j)t{Po) = 0, otherwise, when xq = 0, we have 4't{Po) = ^ffe^o R''"^ ■ In 
any case we have 

MPo) ^ = -F^(Z?20(Po),7^0(Po)) 

and (12.31) is satisfied. 



Assume now touches u from below at Pq. Wc claim that the function can not be touched 
from below by a function at a point (0, to)- If this were the case we would have 

Cht-^A{\xlt)^ -Chto^R^ ^ Bh{x,t) - Bh{Po) ^ t) - 0(0, to) 

> (f>t{x, t){t - to) + -iD^(l){Po) x)-x + o(|a:|2) + o{t) as (x, t) ^ Pq. 

Taking t ~ to and x — Xe where A > and e is any unit vector in M*^ wc obtain 

c,,to* (^(A,to)^ -i?*^) ^ y(-D'0(Po)e) -e + olA^) as(x,t)->Po. 

Since 

^. A(A,to)'^-i - R~ 

lim ^ = — oo 

AiO A^ 

this leads to a contradiction and the claim is proved. 

It follows that Pq must satisfy A(|a;o|,to) = 0. At points where this occurs we have (f>t{Po) = 
dtBhiPo) = and D^{Po) = DBh{Po) = 0, hence 0t(Po) + Fh{DmPo), D(t>{Po)) = and ^ is 
again satisfied. □ 

Remark. Observe that for h G (1,2] the only points where Bh is not twice diffcrcntiable are of the 
form (0,t), hence our proof can be simplified in these cases. 

3.2 Separation of variables 

Let us now look for solutions of the form 

u{x,t)^T{t)Xir), 
where r = \x\. From equation (11.11) we must have 

r \X'\''-^X" 



\Th-i\T X 



-m 



for some to G M. Using the second symmetry in (|2.4p . we can rescale this equation to choose the 
value of TO that better suits us. It is convenient to take to = ±— !-r. 

h—l 

3.2.1 Case to = 

We can immediately integrate T and obtain for an arbitrary to G 

m = - — 

(t-to)^ 
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To integrate X we use the change of variable 

r(s) = K / (sin((T))" da, 



where a = and = 2a (a similar change of variable is performed in jAJKj in the case 

h = 3). Since r'[s) > for s g (0,7r), this transformation is invcrtiblc from [0,7r] to [0,i?], where 



R ^ K (sin((T))" d<T. 
Jo 



Hence, its inverse function s : [0, R] — !> [0, tt] is well defined, strictly increasing, s(0) = 0, s{R) — w 
and s'(r) = k,-^ (sin (.s(r)))"". Let us define : [0,R] R by 

X^{r) := cos(s(r)). 

We have 

Xiir) = --sini-"(,s(r)), X'Jir) = -i^ sin-2"(s(r)) cos(5(r)). 
Since a G (0, 1), the function X^, is [0, R] and C^(0, i?). We can extend X^, to Rq by reflection, 

(3.2) X(r) 



X^{r - 2fci?) if r e [2fci?, {2k + l)R) for fc £ No, 
X^{2kR - r) if r e [(2fc - 2fci?) for k G N. 



Then X is 2i?-periodic function as is the function \X'\^ ^X' . Furthermore, the latter satisfies 

1 - a 



[\X\r)\^''X'{r)]' = -^X{r) = -mX{r) 



even at the points where X is not twice differcntiablc, that is, the points of the form kR, with 
k e No. 

Proposition 3.2. For any to gM. and any tq > 0, the function 
(3.3) 5(x,t;ro,io) = -^^^^, 

{t-to)~ 

h + l 

where Xr{x) = ~ "'^('^Ir^)' with X given by p.2p . is a viscosity solution of (|l.ip m 

M'' X (tQ,(x). In particular, the restriction of S to Qro.to(O) ^'ro(O) ^ (^OiOo) ^ positive 
solution of the homogeneous Dirichlet problem in this set. 



Proof. Using the x-symmctry in (|2.4p . it is enough to consider ro ~ R/2. From the properties 
of X deduced above, it is clear that 5 is a classical, and consequently a viscosity solution of 
(|l.ip off the set E := {{x,t) \ \x\ = kR, k S No}. Moreover, S is in the whole space with 

St{x,t) = -{t-ta)-^^X{\x\) and DS{x,t) - (t ^ t^)-^^ X' {\x\)^y 

We need to check that (|2.3p is satisfied on E. To this effect, let = (x*,t*) G E, > toj 
and suppose (p € C^'^ touches S" from above at P*. As in the proof of Proposition 13. 1[ since 
X"{r) = +00 at the points of the form r = {2k + 1)P, k G No, we see that |a;*| can not be of this 
form, hence — 2kR for some k G No. Now we have D(f>{P^,) = DS{P^,) = and 

c^t{P^) + F{D^p{P,)) = -(<, - to)"^ < 0, 
which is (|2.3p . The verification of the other half of Definition [^T] is similar. □ 
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3.2.2 Case m = --^^ 

In this case we obtain a blow-up solution with 

Tit) = (to-t)-^, 

for t < to^ and 

h + l 

X{r) = ChT'^-^, 

with Ch as in (j3.ip . In fact, we will need a slightly more general solution. 
Proposition 3.3. For any ro,to G K, the function 

V[x,t;rQ,to) = — 

(ta-t) — 

is a viscosity solution of in M'' x (— oo,to)- 

The proof is similar to the proof of the previous proposition, having as the difficult points the 
sphere = tq} when tq > 0. 

3.3 Traveling waves 

Let us finally look for solutions of the form 

u{x,t) = f{xi - ct). 

The function must satisfy 

l/'l'-V7" = c. 

Taking the constants of integration equal to zero (we can recover them using the symmetries in 
(|2.4p ). and ignoring the stationary solutions (c = 0), we integrate twice to obtain 

\c\h 

Therefore, we obtain 

u{x,t)^^{ch^cx^)^\ 
\c\ 

with dh = — — valid for ^ ^ t. Since this function is up to the hyperplane xi = cf, with 
Ut and Du vanishing there, we can proceed as in the two previous propositions to get the following 
one. 

Proposition 3.4. For any unitary vector v and c G K, the function 

dh ^ 
T(x, t] i/,c) — — — cx ■ lA 

is a viscosity solution of in R'* x R. 
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4 Existence 



The main purpose of this section is to prove a solution to (|l.ip - (|1.2p exists. However, as we mention 
in the introduction, and without much additional cost, it will useful to allow a zeroth-order term 
as in (jl.3|) . With e ^ and (5 > 0, we start with the approximation 



(4.1) 

where the d x d matrix A'^'^{p) is given by 



Ut = C-'^{u) A^^^iDu) : D^u + H{u) in Q 
u{P) = g^'*'(P), on dpQ. 



A^'^ip) = el+ (IpP + (52)(''-3)/2p p. 

The proofs of the estimates for the approximate problem (|4.ip are very similar to the corresponding 
proofs of existence in |PVj or jJKj . For clarity we will present them schematically (see also |ASj . 
where a slightly different approximation is used). 



4.1 Lipschitz estimate in time 

We start with Lipschitz regularity in t and then prove the regularity in x. 

Theorem 4.1. Suppose g G C^{dpQ) and u = u^'^ is a smooth solution of (j4.ip . Then there 
exists Ki > depending only on \\D'^g\\oo, \\Dg\\oo, \\g\\oo and ||5t||oo such that 

\u(x,t) -5(2;, 0)1 ^ Kit 

for any {x,t) € Q. If g is only continuous in x and bounded in t then the modulus of continuity of 
u on fl X [0,t*] (for small t^) can he estimated in terms of \\g\\oo and the modulus of continuity of 
go :=5|g^Q inx. 

Proof. Assume for the moment g ^ and H = 0. Let A > and v^{x,t) = go{x) ± Xt. Then, 
if we choose A large enough, u+, respectively v~ , becomes a super-, respectively subsolution of 
(|4.ip which lies above, respectively below u on dpQ. Therefore, by the classical comparison (see 
for example |LSU| ). 

\u{x,t) - go{x)\ < Kot 

for t G [0, T] where Kq is a constant depending only on the stated norms of g. 

If iJ ^ 0, then is still a super/subsolution of (|4.ip in i7 x (0,r*), above/below g on dpQ, if 
we further choose A to satisfy A > 2M|j(7||oo and T* < (2M)~^. Therefore the estimate is valid on 
this time interval. Then we can iterate the argument until we cover the whole interval [0,r]. 

Assume now g is only continuous in x and bounded in t and let loq be the modulus of continuity 
of go. Let us fix a point xq £ fl and < p < min(dist(a;o, dH,), 2y/\\g\\ao)- Let us also define 

g^{x,t) = go{xo) ± wo(p) ± \x - xop. 

P 

It is easy to see that g~ ^ ^ g+ on F and thus, again from the comparison principle, u~ ^ f ^ u"*", 
where is the solution of (|4.ip with initial and boundary condition g^ . Since g^ are in (M.'^ x M) , 
we can use the above estimate to conclude that 

|u=^(a;o,t) - g^(x-o)| ^ K^t 
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where depends on ||.9[|oo and p. Therefore, 

3 

\u{xo,t) ~ goixo)\ < 2Kf^t + -uJo{p). 
This inequahty concludes the proof. □ 
Using (j4.ip we obtain the fuh Lipschitz estimate in time. 

Theorem 4.2. If u is a solution of (|4.ip in Q and g G C^{Q), then there exists A'2 > depending 
only on \\D'^g\\oo, \\Dg\\^, \\g\\oo and \\gt\\oo such that 

\u{x, t) - u{x, s)\ ^ K2\t - s\ 

for every x t, s Cz (0,r). If g is merely continuous, we can estimate the modulus of continuity 
of u on Q in terms of ||.9o||oo and the modulus of continuity of ga. 

Proof. Taking t > and 

t) := u{x, t + t), 
using the Theorem 14.11 it is immediate to get 

\u{x,t) ~ uix,t)\ ^ K2t 

in fl X [0, r — t]. The case when g is only continuous is done as in the previous proof. □ 
4.2 Holder continuity in space 

Theorem 4.3. Let u be the solution of (|4.1|) with g G C^((5) nLip(Q). There exist a £ (0, 1) and 
Kt, > 1, depending only on, ||.g||cx3, ll-DjfUoo, HfftHoo and a, such that for every s and 6 sufficiently 
small and for every Pq ~ (.tq, to) £ T and x G f2 with \x — Xq\ ^ \ we have 

|M(x,io) - g(a;o,io)| ^ K'i\x - xoT- 

Proof. Let us define 

v^{x,t) = g{xo,ta) + K^\x - .to|" + A (to - , 

where K* 1 and A > are constants which we will choose in such a way as to make v'^ a 
supersolution lying above u on the appropriate domain. Let us take x G rin-Bi(xo) and t G (0,to)- 
An easy computation, using the fact that jx — xo| ^ 1, yields 

where 

Since jcc — xo| ^ 1, we need only consider z ^ aK^,. If h ^ 3, then f{z) ^ and when 1 < h < 3 
we have f{z) > 2^'^~^'>/^z'^~^ ^ z''"^ provided z > 5. In any case it holds 

- > -A + f '„"J p-^At-^ - - A/(i|5||oo + + TA) > 
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whenever e,S < 1 < aK^, 

aK^ > (^^Y~ir) " ' '^^^ ^^ll.9lloo + (1 + MT)X < (1 ~ a)a^Kl^ - {e{d - 1) + M)aK^. 

We want to have ij+ > u on Q* = (il n Si (xq)) x (to — t*,to)- where we take := min{l,to}- Let 
P = {x, t) G dpQ* . Let us first assume P is on the lateral boundary of Q* . Since u ~ g on dQ and 

\x - Xo\ < 1, 

uiP) ^ uiPo) + \\Dg\\^\x ~ xo\ + ||5t||oo(<o - t) 
^ g{Po) + K,\x~ xol" + \{to - t) = v+{P), 

provided A'* ^ ||D(7||oo and A ^ ||.gt|ioo- If, on the other hand, x e n dBi{xo), then, using 
comparison, 

u{P) ^ II5II00 s$ u{Po) + lu + A(to -t)^ v+iP), 

provided > ||ff||oo- 

We consider now the case when P is on the bottom of the cylinder Q*. Let us first assume 
X E fin Bi(xq) and t — to — I. In this case, again using comparison, we get 

uiP) ^ Moo ^ u{Po) + R\\x - xo\" + \ = v+iP), 

as long as A ^ llff lloo- Finally, when to < 1, and hence Q* ~ {flD Bi) x (0, to), we have that u = g 
on the bottom, therefore 

u{P) = u{x,0) = gix,0) < g{xo,tQ) + \\Dg\\oo\x - xa\ + \\gt\\ooto 
^g{Po) + K^Ax-xor + Xto^v+{P), 

provided, once again, ^ ||D(7||tx3 and A ^ ||5t||oo- 

To summarize, we have ^ m on dpQ* , and hence, by comparison, f ^ ^ u in Q* , if e,(5 < 1, 

K^a ^ maxjl, HA^Iloo, llffll 

^llslloo + (1 + MT)\ < (1 - a)a^K'l - {e{d - 1) + and 
A ^ max{||gt||oo, llfflloo} • 

Since /i > 1 we can choose the constants K^, and A to satisfy these inequalities. We obtain 

u(x,to) - 9ixo,to) v+{x,to) - g{Po) = K^\x - xo\". 

Using the barrier :~ g{Po) ^ K^,\x — xqI" + X{t — tp) we get the reverse inequality, 

u{x,to) - gixo,to) ^ -Kt\x - xo\". 

□ 

We can extend the estimate to the interior of the domain. We will use the following notation for 
convenience. For z S M'', we define il.z = z + n = {x + z\ xGft} and tor r > 0, ilr = {x (£ \ 
dist(a;, dV,) ^ r}. Given x, ?/ S M'^ we define the closed segment [a;, y] := {6y+{l — 9)x | ^ 6* ^ 1}. 
The semi-open and open segments [x,y), {x,y] and {x,y) are defined analogously. 



e{d-l) 
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Theorem 4.4. The conclusion of Theorem \4.!^ is valid in the interior of Q, that is, there exists 
K4, depending only on \\g\\oo, ||^.9||oo cind ||5t||oo; such that for e and S sufficiently small and for 
every x^y € fl 

\u{x,t) - u{y,t)\ < Ki\x - 2/1". 

Proof. To simplify notation, in this proof wc omit the dependence of u in t, since time does not 
play any role. Take a vector z G i?i(0) and let F = $7 n $7^. Let us define Uz{x) := u{x — z). From 
Theorem [43] we have that \u{x) — Uz{x)\ ^ i^slzj" on dV (x S dV implies that x G d^l or x — z € 
dn). Hence, using the comparison principle we have that Uz{x) — K^lzl" ^ u{x) ^ Uz{x) + /Cslzj" 
for X £ V. This means that whenever x,y G ^x-y H f2 or x, y G ^y-x H fi, with \x — y\ ^ pa, we 
have \u{x) — u{y)\ ^ K^lx — ?/|". In particular, the same is true whenever x,y € Q,\x-y\- 

When \x — y\ > 1, using the comparison principle we obtain the conclusion of the theorem taking 
A'4 ~ 2||g||oo- Let us therefore assume that |x— y| ^ 1 and x~y ^ Let us first further assume 

that [x,y] C fl. In this case we can take the two segments [x, and [w,y], where w = {x + y)/2 
is the midpoint of [x, y], let z = y — w and note that w,y G flz r\fl and x,w & ^-z H ri. Hence, 
from the first step of this proof, we have 

\u{x) — u{y)\ ^ \u{x) — u(w)\ + \u(w) — u{y)\ ^ K^{\x — wj" + jui — y\°') 
2i""A'3|a;-y|". 

If the segment [x^y] is not completely in VL, then we can certainly find 101,1(72 £ dVL n [x,y\ 
(not necessarily different) such that [x,wi) G and {w2,y\ G fi. In this case we can apply 
Theorem 14 . 31 directly to get \u{x) —u{wi)\ ^ i^Tajx — and \u{w2) —u{y)\ < K3\w2 — Since 
\u{wi) — u{w2)\ ^ ||Dg|||wi — W2I ^ K^lwi — W2\°', we easily get the result with K4 = A^^'^K^. 
This finishes the proof. □ 

4.3 Lipschitz estimate in space (e = 0) 

As in [PVj , we obtain Lipschitz regularity in space when we take e = in the approximation (|4.ip . 

Theorem 4.5. Let g G Lip((5) and suppose u is a viscosity solution of ()4.ip with e = 0. There 
exist a constant K^, depending only on \\g\\oo, ||I?.g||oo and ||5t||oo (independent of S G (0,l)j, such 
that for every Pq = ixo,to) G dfl x (0, T) and x € flCi Bi{xo) we have 

\u{x,to) - g{xa,ta)\ ^ K^lx - xo\. 

Furthermore, if g is only continuous, then the modulus of continuity of u can he estimated in terms 
of \\g\\oo and the modulus of continuity of g. 

Proof. Let A'*, and A be positive constants and define 

w+(a;,t) = g{Po) + L^\x - xo\ - K^\x ~ xq]'^ + \{to -t). 

We will check u+ is a viscosity strict supersolution on Q* {Q* defined as in the proof of Theorem l4.3p 
above g for an appropriate choice of constants L*, and A. Observe that 

vt{x,t) ~ C"'^v+ix,t) = -X + 2KJi\L,-2K,\x-XQ\\)-H{v+{x,t)), 

where, as in the proof of Theorem 14. 3[ f{z) = (z^ + S^y^^^^^^z^ . Recall from that proof that 
f{z) > z^~^ for z ^ 1 and for all h > 1. Hence, if we choose and AT* such that > 3A'* and 
2jj^^+i _ 2MK^ > (1 + MT)X + MII5II00, the above inequality shows that is a supersolution. 
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We need to further choose the constants so that w+ ^ / on the parabohc boundary of Q* . Let 
P = {x,t) be a pomt in F* = dpQ* . If a: e dfl, as before 

uiP) = g{P) ^ .g(Po) + \\Dg\\oo\x -xo\ + \\gt\\oc{to - t) 
sC g{Pa) + [L, - l)\x -xo\ + \{to - t) < v+{P), 

provided ^ HDpHoo + 1, K^, ^ 1 and A > ||gt||oo- If a; G H dBi{xo), then 

u{P) ^ biloo ^ g{Po) + {L, - l)\x - xol + A(<o - t) < v+{P), 

provided L* ^ llslloo + 1 and yet again ^ 1. 

When P = {x,t) is on the bottom of the cyUnder Q*, as before we consider two cases. When 
to ^ 1, t = to — 1 and 

uiP) ^ llfflloo <.g(Po) + A<f+(P) 

as long as A ^ ||.9t||oo and ^ max{l, A'*}. On the other hand, if to < 1, t = and hence 

u{p) = g{P) ^ ||.g|U g{Pa) + A < v+{P) 
under the exact same conditions as for the previous formula. 

Therefore, we have u ^ w+ on F*, and thus by comparison on Q*, as long as we take 

A^||.gt||oo, if* ;5max|l,v/A72} and ^ max {2, + 1, !|.9||oo + 1, 3/C} • 

Using once more the comparison principle, we have that for x E^C] Bi(xo), 

u{x,to) s$ v'^{x,to) ^ g{Po) + L^\x - xo\. 

Using instead the barriers 

v~{x,t) = g{Po) - L^\x- xo\ + K*\x - xo\^ + A(t - to) 
we obtain the reverse inequality and, as a consequence, the Lipschitz estimate. 

Let us finally merely assume that g is continuous and let Wg(cr) be a modulus of continuity at Po- 
More specifically, let ujg be a continuous, decreasing function in a such that \g{P) ~ <7(-Po)| ^ ^g{(^) 
whenever max{|a; — xo|, |t — to|} ^ cr. Let a g (0,to) and define the smooth functions 

g^{x,t) := .g(2;o,0) ± Wg(o-) ± 5— a; - xo\ ± \t ~ to . 

If max{|x — xqI, |t — toll ^ cr, then 

g-{P)^g{Po)-Ug{a)i^g{P)^g{Po)+ujg{a)^g+{P), 
and if max{|x — xol, |t — to|} ^ ct then 

g^{P) ^ -II.9II00 g{P) ^ \\g\\oc. ^ g+{P). 

Therefore, if are the solutions of (j4.ip with e = and initial data g*, by comparison ^ u ^ 
on Q. Since are smooth we can apply the first part of the theorem to deduce that 

\u^{x,to)~g^m\^K+\x~xol 

where K'^ depends on H^Hoo and a. From these inequalities we get 

3 

\u{x,to) ~ g{Po)\ < 2K+\x~xo\ + -^^ojgicj). 
This finishes the proof. □ 
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Our final estimate is the interior Lipschitz estimate. 



Theorem 4.6. Let g and u he as in Theorem \4.5\ For every x,y £ D, and t G (0,r) 



ti(a;, t) 



u{y,t)\ < A'six - y|, 



where is the constant given in that theorem. If g is only continuous, then the modulus of 
continuity of x i— >■ u{x,t) can be estimated in terms of \\g\\oo cind the modulus of continuity of g in 

X. 

Proof. The proof is similar to the proof of Theorem 14.41 but in this case it is easy to get the 
optimal Lipschitz constant. Once again we omit the time dependence of u. Take z 6 K'' such that 
\z\ ps. Define V = Cl flz and let Uz{x) := u{x — z). From previous theorem we know that 
\u{x)—Uz{x)\ K^\z\ ondV. Using comparison, we have that Mz(a;)—iir5|z| u{x) ^ Uz{x)+Kq\z\ 
in V. Therefore, \u{x) — u{y)\ ^ Kq\x — y\ H x,y G ^x-y, and in particular the same is true if 
x,y £ n^x-y\, Or = {x e I dist(a;, 917) > r}. 

Suppose now x — y ^ ft^^-yi and let us first assume that the whole segment [x, = {z G M'' | z ~ 
9y + {1 — 0)x, ^ 6* ^ 1} is in VI. Let us assume without loss of generality that p = dist(a;, dVl) ^ 
dist(j/, We can find points Xi, Q ^ i ^ n such that x = xq, Xn = y, Xi € [x^-i, x^+i] 
(1 ^ I ^ n — 1), and pi = \xi — Xi-i\ ^ p (1 ^ i ^ n). Noting that Xi^Xi-i G Vlxi-xi^i we can use 
the previous step to conclude that \u{xi) — u{xi^i)\ ^ K^\xi — x^-ij, and hence 



If, on the other hand [x, y] ^ i7, then we can find points xi,X2 G dO. n [x, y] such that [a;, xi) dO. 
and [x2,y\\{x2\ C O. We can further choose wi G [x, xi], with |wi — xi | ^ 1 and W2 G [x2, y], with 
\w2 — X2I ^ 1. Then we apply the above to obtain \u(x) — u(wi)\ ^ K^\x — wi\, \u{w2) — w(y)| ^ 
\w2 — y\, while from the previous theorem, \u{wi) — u{xi)\ ^ Kc,\wi — Xi\. Putting all these 
inequalities together gives the Lipschitz estimate for this last case. 

The proof of the statement with the modulus of continuity follows as in the proof of Theorem l4.4l 



We finally prove Theorem 11.11 Existence is proved by piecing out the results in Theorems 14. 2[ 
14.41 and 14. 6i and using the standard compactness arguments, as is done in |PV1 Theorem 1.1]. 
Uniqueness follows directly from the comparison principle, Theorem 12.21 

Proof of Theorem ] l.U The proof when is bounded is similar to the proof of jPV[ Theorem 1.1] 
and we just sketch it here. Assume first g G C^{Q) H Lip((5). The comparison principle and 
Theorems 14.21 and 14.41 imply that the family of functions {u"^'^} is uniformly bounded and equicon- 
tinuous, therefore, for some sequence Sk — >■ 0, u'^'"^ , which, by a standard argument of 

viscosity solutions, is a solution of (|4.ip with eps ~ 0. Then, using Theorems 14.21 and 14.61 we can 
find a sequence f^'' — !> u. The stability arguments for viscosity solutions work here to show that u 
is a viscosity subsolution of (|l.ll) . 

If n is not bounded we define nR:=nn 5^(0), Qr := Qr x (0, T] and gR-.TR-^ M, = dpQR, 
by gR{x,t) := if |x| = R, gR{x,t) := XR{x)g{x,t) for (x,t) G F n Qr, where xr{x) = xi^/R) 
and X 6 C^iR'^) satisfies x(x) = 1 if |x| ^ 1/2, x(x) = if \x\ ^ 1. In Qr there exists a unique 
viscosity solution with initial/boundary data gR. From the assumptions on gR, the estimates for 
ur and the stability of the viscosity solutions we can let i? — )• 00 to obtain the result. □ 



n 




□ 
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5 Asymptotic behaviour in the whole space 



We consider in this section Cauchy problem (jl.4[) . We first obtain a decay rate for the solutions of 
this problem and then prove the asymptotic convergence. 

5.1 Decay rate 

Using the similarity solutions and the comparison principle we readily obtain the following esti- 
mates. 

Theorem 5.1. Let u he the unique viscosity solution of (|1.4p with uq G Cc(IR'^) not identically 
zero, Uq ^ 0. Then, there exist positive constants c and C depending only on Uq such that for t > 

(5.1) c(l + t)"* ^ max|M(a;,<)| s: C(l + t)-^. 

Moreover, the support of u expands continuously at a rate of the order of t^ . 

Proof. We can find similarity solutions Bi and B2 of the form (|3.ip such that 

±Bi(x,0) ^ uo{x) s: ±B2{x,0), 

which immediately gives (|5.ip . Indeed, if Xq is a point where uq(xo) ^ 0, then, with r sufficiently 
small we can take -Bi(a;, t) — i3r,/i(a;—xo, t+1) and with i? sufficiently large -82(2:^, = Bii,h{x,t+1). 

To prove that the support expands continuously, take a point on the boundary of the support 
of u and the take a function of the form Bi, as above, with Bi ^ u, with center sufficiently close 
to P*. Then, the comparison principle implies that the support cannot jump "inward" toward the 
support of Bi. In the other direction, we can find a function V2 of the form given in Proposition l3.3l 
such that w ^ V2 for a short time. Again the comparison principle implies that the support can not 
jump "outward" beyond the support of V2. The rate of expansion of the support of u is controlled 
by the rates of expansion for Bi and B2 above, that is, it has to be t^h. □ 

5.2 Asymptotic behaviour. Proof of Theorem 11.21 

Step 1. The idea of the proof is the same as in the proof of |PV1 Theorem 1.4]. First we let Bi 
and B2 sandwich u as in the proof of Theorem 15. II Then consider the family of rescaled solutions 
u^, Bi and B2, where for a function f{x,t) we define f^{x,t) by 

f^ix,t)^X^^f{\^x,Xt). 

Since Bi and B2 are invariant under this transformation, on any compact time interval [ti, t2] with 
< ti < t2 < 00 the family is continuous, uniformly bounded, and supported on a uniform ball 

Step 2. Now we use Aleksandrov's principle, as explained for instance in |CVW| . to show that for 
a solution u(x, t) with initial data uq{x) ^ supported in the ball Pfl(O) we have for alH ^ and 
all r > P 

inf u(x,t) = max u{x,t) 

\x\—r \x\—r~^2R 
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Note that in doing this, we need to use the traveling wave solutions from Proposition 13.41 to show 
that the solutions are almost radial. If this is applied to the rescaled solutions, we get for all 

|X| >i?A =i?A-l/(2'') 

inf u^{x,t) = max u^{x,t). 

\x\—r \x\—r~\-2B.x 

Step 3. We now fix t = 1, A very large, so that R\ ^Si e is very small, and define 
ui{r) := inf u^{x,l), U2{r) := maxu^{x,l). 

\x\—r \x\—r 

We easily verify that ui (r) , U2 (r) are nonnegative and radially symmetric functions, both supported 
in the same ball _B/f^(0), they are nonincreasing as functions of r for r ^ e, and we also have 

U2{r) > ui{r) > ?ii(r + e) 

for all r ^ e. It is then easy to verify that the 1-d mass of U2{r) — ui{r) is less than Ce. 

Step 4. If ui{r,t) and U2{r,t) arc the corresponding radial solutions of the problem with initial 
data at t = 1 given by ui{r) and U2{r)^ respectively, we have for all t ^ 1 

ui(r,t) u^{x,t) ^ U2{r,t) 

As with the convergence result for the 1-d PME, cf. |Vl Theorem 18.1], the result follows. 

6 Asymptotic behaviour on bounded domains 

In this section we analyze the asymptotic behaviour of the solution of the homogeneous Dirichlet 
problem (jl.Sp . As in the previous section, we first obtain a decay rate for the solutions of the 
problem and then prove the asymptotic convergence. 

6.1 Decay rate 

We start by proving that the decay rate for the solutions of (jl.5|) is t^T^ , the same as the decay 
rate of the friendly giants in p.3p . Of course, this is not a coincidence. 

Theorem 6.1. Let uq satisfy the condition of Theorem ] 1.31 let u be the unique solution of (jl.Sp 
and assume Xo is a point where uq(xq) > 0. There exist positive constants to, ri and r2 such that, 
with Xr as in (|3.3p . there holds for every t > 0, 

Xr-^{x — Xq) !^ {t + to)''^ u(x,t), for X € Brj^{xo) and 

u{x,t) ^ Xr2{x — xo), forx£Q. 

Proof. It it straightforward, under the conditions of the theorem, to find ri and to such that 
S{x — Xq, 0; ri, —to) ^ uq{x) on the ball (a;o) and ^2 such that ft C (^^o)- The result is then 
immediate from the comparison principle. □ 



17 



6.2 Friendly giant in Q. Proof of Theorem 11.31 



Let us now show there exists a friendly giant in 51. We will show that the asymptotie profile for 
arbitrary (nonnegative) initial conditions is the profile of this friendly giant. 

According to (12.41) . we can reseale a solution u of p.5p by 

u\{x,t) :~ A*"-! u{x, At) 

and still obtain a solution of the same problem but with initial condition uox{x) = A''^uo(a;). For 
A < 1, by comparison, we see that u\(x,t) 5C u{x,t). According to [BCl Theorem 2.3], we deduce 
that 



(6.1) 



t + t) — u{x, t) ^ — 



1 - 



t 



t + T 



u[x, t) 



for (.T, t) G Q and r > such that t + t < T. Let us now consider the following rescaling of u, 

v{x, s) {h - l)T^e'u (x, e^'-^)'^) . 

It is easy to see that w is a viscosity solution of 

' Vs — ^'^v = V inr2x(l,oo], 
v{x,0) ^ {h - l)'^u{x,l) {orxeil, 
v{x, s) = for x <E dQ, s > 1. 

The estimates for u in Theorem 16. II and in (j6.ip imply the following estimates for v: 

(6.2) v{x,t)s^M, v{x,s + h) -v{x,s) ^0. 

Using our similarity solutions p.ip . the proof of jLS| Lemma 3.2] adapts step by step, mutatis 
mutandis, and we have that v eventually becomes positive on any compact subset of 57. 

Lemma 6.2 ( [LSI Lemma 3.2]). For any compact set K C there exist sk CLud mx such that 

v{x, s) ^ ruK on K x [sk, oo). 

Using this and the estimates (|6.2p . we deduce that there exists a lower semicontinuous function 
Go : 17 — )- M such that Gn{x) = on dH. and 

lim v{x, s) = Gfiix). 

s— f oo 

Since w is a viscosity solution of (|1.3p . from Theorem ll.il we see that we can control the modulus 
of continuity of v and so Gn must be continuous. In terms of u, this means 

lim tT^u{x,t) = Fn{x) := {h - 1)-t^ Gn{x). 
Theorem 6.3. The function Gn 'is a positive viscosity solution of the eigenvalue problem 

-A^Go = Go m n, 
Go(a:)=0, forxedn. 

The proof is a simplification of [PVi Theorem 7.3]. With this we finish the proof of Theorem ll.3l 
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7 Comments and extensions 



• The restriction of nonncgativity in the study of asymptotic behaviour can be somewhat weakened 
by observation that the equation is invariant under constant displacement of the u- variable. Hence, 
we can assume that the initial data arc compactly suported perturbations of the level it = c, c 
constant. On the other hand, the restriction to compactly supported data can be weakened, but 
we do not know how to find a reasonably wider class, much less an optimal class. 

• The solutions of the form u[x,t) = t^'^^ Xr{x) described in p.3p are interesting examples of 
radial solutions which decay in time but not in |a;|; they oscillate radially in a sine-like fashion. 
This behaviour is distinct from the several-dimensional heat equation, where it is known that such 
solutions decay like a Bessel function as goes to infinity. The restriction of our functions to 
appropriately chosen balls also gives an example of an asymptotic profile which changes sign as 
many times as desired. This naturally leads to the question of classifying asymptotic profiles for 
sign-changing initial data. We do not know if it is possible to obtain nonradial stable profiles on 
a ball. 

• In both settings, Dirichlet problem in a bounded domain with zero boundary data and Cauchy 
Problem in the whole data, we generate a semigroup enjoying the comparison principle. Since the 
operator is h homogeneous with h ^ 1^ the homogeneity estimate of Benilan-Crandall [BC| applies 
and we have for all nonnegative solutions 

in the sense of distributions. This estimate could be important in deeper studies. 

• The "zeroth-order" large time behaviour for the Dirichlet problem with nonhomogeneous bound- 
ary conditions is considered in |AJKj . We do not consider this case here, but the techniques used 
there should work for our class of operators. 

• We are not exploring the cases h = \ which is a well-known equation, or the cases h < \ which 
must have novel properties. 
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